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f^ i Abstract. The phase diagram of QCD, as a function of temperature T and quark 

^N I chemical potential fi, may contain a critical point {ij,e,Te) whose non-perturbative 

nature makes it a natural object of lattice studies. However, the sign problem prevents 
the application of standard Monte Carlo techniques at non-zero baryon density. We 

iy-\ I have been pursuing an approach free of the sign problem, where the chemical potential 

is taken as imaginary and the results are Taylor-expanded in /z/T about /i = 0, then 
analytically continued to real fi. 

Within this approach we have determined the sensitivity of the critical chemical 
(^ . potential ^ie to the quark mass, d{^E)'^ /dmq\fj^^=o. Our study indicates that the 

D I critical point moves to smaller chemical potential as the quark mass increases. This 

f^ , ' finding, contrary to common wisdom, implies that the deconfinement crossover, which 

takes place in QCD at ^ = when the temperature is raised, will remain a crossover 
in the ^-region where our Taylor expansion can be trusted. If this result, obtained on 
a coarse lattice, is confirmed by simulations on finer lattices now in progress, then we 

\^ ', predict that no chiral critical point will be found for hb ^ 500 MeV, unless the phase 

OO ' diagram contains additional transitions. 

o 

00 . 

(^ , 1. Motivation 

QCD has a rich phase diagram as a function of quark chemical potential fi and 
^ , tem-perature T, with at least 3 regimes: confining (low fi,T), quark-gluon plasma (high 

T) and color superconducting (high /i) . Lattice simulations have provided clear evidence 
that, for /i=0, the temperature-driven "transition" between the first 2 regimes is actually 
a crossover [1]. Then, the commonly expected (/i,T) phase diagram, depicted Fig. 1 
(left), contains a critical point caused by the /x = crossover turning into a first-order 
transition. This critical point is the object of both experimental search and theoretical 
lattice investigation. If located at small enough chemical potential, it could be discovered 
in high-energy, not-so-heavy ion collisions at RHIC or LHC. In fact, it has been predicted 
to lie at (/is, Ts) = (120(13), 162(2)) MeV in a celebrated lattice study [2]. 

However, the numerical method of Ref. [2], which consists of reweighting results 
obtained from /i = simulations to /i 7^ 0, is known to fail and give wrong results without 
warning (see, e.g., [3]) when the overlap between the /i = Monte Carlo ensemble and 
the target /i 7^ ensemble becomes insufficient. Therefore, the abrupt change observed 
in [2] as a function of fi (see Fig. 1 there) might be caused by an abrupt breakdown of 
the numerical approach. In addition the needed reweighting factors to the fi ^ target 
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Figure 1. Conventional (left) and simplest possible unconventional (right) phase 
diagram of QCD as a fmiction of chemical potential /i and temperature T. Dotted 
lines denote crossovers, solid lines first-order phase transitions. The QCD critical 
point (end of first-order line), marked by a 'v' on the left, is absent on the right. 



ensemble are not always positive, because of the notorious "sign problem" affecting the 
/i 7^ determinant. This makes the statistical errors more difficult to control. Finally, 
the lattices used in [2] have Nt = 4 time-slices, corresponding to a very coarse lattice 
spacing a ~ 0.3 fm. As acknowledged in [2] already, the results may change considerably 
after extrapolation to the continuum limit a — > 0. Therefore, the landmark result of [2] 
should not be considered the final word on the QCD critical point. A more cautious 
crosscheck like the one we present here is warranted. 

On the theoretical side, the expectation of a critical point as in Fig. 1 (left) can be 
traced back to two assumptions: (z) QCD with Nf = 2 massless flavours, at fi=0, has 
a second-order temperature-driven phase transition in the 0(4) universality class [1]. 
This implies, for this theory, the existence of a tricritical point at {fi, T) , and for the 
Nf=2 + 1 theory that of another tricritical point at (/x=0,T*,m„,i = 0,m*). (m) These 
two tricritical points are analytically connected by a line in the (/i, T, mu^d = 0, ""^s) space. 
However, assumption (i) depends on the strength of the axial Ua{^) symmetry breaking 
at Tc. And assumption (ii) has no other basis than simplicity. This motivates us to take 
a careful look at the phase diagram of Nf=2 + 1 QCD, generalized to arbitrary quark 
masses mu=md=mu,d and mg. We present our /i = 0, then fJ^^O findings below. 



2. Approach 

To protect ourselves against the pitfalls of reweighting to /x 7^ 0, we consider the effect 
of an imaginary chemical potential /i = i^j [71 [6]. Then the fermion determinant is 
positive and no overlap problem occurs because we do not reweight. The drawback is 
that, in order to analytically continue an observable, its /i-dependence is fitted to low- 
order Taylor series about yU = 0, thus introducing a truncation error. As a safeguard, we 
also calculate the coefficients without truncations [9]. Note that these derivatives can 
be expressed, and measured, as non-local observables in the fi = ensemble. This is the 
strategy pursued successfully in [8] for the pressure. Our approach is computationally 
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Figure 2. Order of the tcmpcraturc-drivcn phase transition (at fi — O), as a function 
of the Hght quark mass iJiu^d and the strange quark mass ttIs- On the left, theoretical 
expectations [5]. On the right, iV^ = 4 lattice results for light quarks, from [6]. The 
straight line corresponds to Nf = 3 (m^ — ra^.d)- The curve is the expected scaling 
behaviour from a tricritical point at m^^d — 0, ttIs ^ 500 MeV. The physical point, 
marked by a x , is in the crossover region. 



more economical, and we extend it to the critical line. 

The first step is to determine, at /i = 0, the line in the {mu^d, i^s) plane for which the 
temperature-driven phase transition is second order. For lighter quarks, the transition 
becomes first order as in the chiral limit; for heavier quarks, lattice simulations show a 
crossover. Theoretical expectations are displayed in Fig. 2 (left). Our determination of 
the critical line on a.n Nt = 4 (a ~ 0.3 fm) lattice is shown Fig. 2 (right), and matches 
expectations, including the possible scaling behaviour in the vicinity of a tricritical point 
at (wT-u ,^ = 0, m*"'^ ~ 500 MeV). Also, the physical point (with real- world it and K meson 
masses) lies, as expected, in the crossover region. 

We now choose a point on the critical line just determined, and monitor the effect 
of a small fi. This exercise has been performed, with similar results, at the two points 
marked by arrows in Fig. 2 (right), with higher accuracy for the Nf = 3 case. If 
conventional wisdom is right, then the surface spanned by the critical line as /i is turned 
on will bend towards the crossover region, so as to produce a critical point at some 
small n value for physical quark masses. This is depicted Fig. 3 (left). On the contrary, 
an opposite curvature will indicate the absence of a critical point at small /i, where our 
truncated Taylor expansion is a good approximation. This is illustrated Fig. 3 (right). 
The issue is thus to determine whether the /i = second order transition becomes first 
order or crossover for a small /i. We have now confirmed our first study [Sj by a second 
one using higher statistics and an improved numerical method [9]. We measure the 
change in the Binder cumulant Bi{X) = {{X - (X))^) / {{X - {X))^ , with X = V^V'. At 
the second order transition, B4 takes the value 1.604 dictated by the 3d Ising universality 
class. Reweighting with a small imaginary /i induces small changes in i?4, which can 
be measured accurately because statistical errors largely cancel between the reweighted 
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Figure 3. Critical surface spanned by the critical line of Fig. 2 as /i is turned on. 
Depending on the curvature of this surface, a QCD critical point is present (left) or 
absent (right). For heavy quarks, the curvature (right) has been determined in [lOj . 
The arrows indicate the effect of a finer lattice: the distance between the physical quark 
masses and the critical surface increases, driving the critical point (left) to larger fi. 



and original ensembles. Our current result for the Nf=3 {niu^d '- 
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supporting the unconventional Fig. 3 (right), both in the leading and subleading term. 
It is now mandatory to repeat this study on a finer lattice, in order to gain some 
control over the continuum limit. Nt = 6 (a ~ 0.2 fm) simulations are in progress. The 
first step, at /i = 0, already reveals an important shift of the critical line towards the 
origin: for the Nf = 3 theory, the pion mass, measured at T = with the critical quark 
mass, decreases from 1.6 Tc to 0.95 Tc [9J. As shown Fig. 3, this considerably increases 
the distance of the critical surface to the physical point, pushing the critical point to 
larger fi (left), or requiring larger higher-order terms of the right sign to bend the critical 
surface "back" (right). Regardless of the sign of the curvature in the continuum limit, 
this trend alone makes a QCD chiral critical point at small fig/T < 1 unlikely. 
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